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We write the averaged Einstein equations in a form suitable for use with Newtonian gauge linear 
perturbation theory and track the size of the modifications to standard Robertson-Walker evolution 
on the largest scales as a function of redshift for both Einstein de-Sitter and ACDM cosmologies. 
In both cases the effective energy density arising from linear perturbations is of the order of 10 -5 
times the matter density, as would be expected, with an effective equation of state iw e g ~ —1/19. 
Employing a modified Halofit code to extend our results to quasilinear scales, we find that, while 
larger, the deviations from Robertson- Walker behaviour remain of the order of 10 -5 . 
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I. INTRODUCTION 

Observations of the cosmic microwave background (CMB) [ESQ, large-scale structure (LSS) (e.g. [1 i]) and 
supernovae type la data [1, 0, @, Q consistently indicate the presence in the universe of significant "dark" components. 
Data from the LSS and from nucleosynthesis indicates that the density of baryonic matter should not exceed about 
5% of the critical density, while data from the CMB suggests that the universe should be flat, implying that around 
95% of the matter-energy content of the universe is unobserved. In combination with numerical studies (e.g. flQj ) 
and observations of the baryonic acoustic oscillations in the LSS [J, |5j this has lead to the "concordance model" in 
which roughly 5% of the universe is in the form of baryonic matter, 20% of the universe is in the form of some "dark 
matter" that interacts only gravitationally, and about 75% is in the form of a "dark energy" fluid with a negative 
equation of state today w « — 1. 

The most popular alternatives to a cosmological constant involve exotic fluids that violate the stron g; en erg y condi- 
tion in the late universe, of which scalar-field models such as quintessence are the most popular (e.g. [Ill [l2l . Il3j for 
pioneering works and [lj, [l5| for recent reviews.) There are also many other modifications of both the matter and 
gravitational sectors. Another approach is to consider the impact of local inhomogeneities on the luminosity distances 
in the local universe (see for example [H, [l?], [H| f° r some recent work on this topic) . The modifications that the local 
inhomogeneities introduce can, it has been suggested, account for many of the features exhibited by dark energy. 

However, the assumptions used to build up our standard model include an implicit averaging procedure, assuming 
that the averaged Einstein tensor is equivalent to the Einstein tensor built from an average metric - and such is not the 
case [ill Hoj]. The Einstein equations are non-linear and local; the correct approach is to average the local equations 
across some domain rather than to assume ab initio "averaged" equations. Studies into averaged cosmologies include 
[2ll . [22I I23I Izil . |25| and the impact of inhomogeneities were applied to account for a dark energy in [26| . In the years 
since, this "backreaction" and the impact of inhomogeneities on observables has been relatively well studied in a 
variety of models and remains an active field HE M, M, M, 11 HE H HE H H3, [H, [H, E3, S3, S3, EE H, SF 



|46L 1471 l48l 149 . |5fJ, [5l|, [52j, [53|; see also [54( for a recent status report. The attraction of this approach is that it can 
recover quantities that - in principle - could act as a dark energy, without the need for exotic matter components 
or modifications to general relativity. Moreover, one can intuitively state that as the universe grows increasingly 
non-linear, the deviations on averaging from a truly homogeneous and isotropic model should increase, potentially 
providing an appealing solution to the coincidence problem. The current universe is significantly inhomogeneous at 
scales of up to 100-300Mpc//i [55|. While we might still expect to recover an FLRW-like metric when averaging across 
scales larger than 100Mpc//i (see [H, H3] for recent studies) there is no automatic guarantee that it will obey the 
familiar Friedmann and Raychaudhuri equations. It has been shown in many studies that structure - both at the 
linear order and for highly inhomogeneous models - can introduce modifications to the large scale effective Friedmann 
equations (though this is still somewhat debated; e.g. [IE Hi], |43|, [58| ) . Even should the effect be relatively minor, 
with current and future experiments a divergence from Robertson- Walker behaviour on the order of 10 -3 could be 



'Electronic address: J.Behrend@thphys.uni-hcidelberg.de 
^Electronic address: I.Brown@thphys.uni-heidelberg.de 
tElcctronic address: G.Robbers@thphys.uni-heidelberg.de 



2 



significant. 

In this paper we calculate numerically as a function of redshift the size of the deviations from standard Robertson- 
Walker behaviour arising from linear and quasilinear perturbations. Such a calculation complements the recent 
studies by [H3|, in which the authors reconstruct the impact of backreaction effects from the observational data, 
and [HH wherein the authors evaluate the size of the effective density of the backreaction as a function of redshift 
in a structured Robertson- Walker model. Much of the current literature concerns exact inhomogeneous rather than 
perturbative models, not least because the impact in a perturbative approach is not expected to be large. However, the 
magnitude and nature of the backreaction and other corrections associated with averaging that arise naturally within 
a standard Einstein de-Sitter or ACDM cosmology remain open questions and this should be addressed. We work in 
conformal Newtonian gauge, in contrast to much of the literature which concerns itself with synchronous gauge, for 
two main reasons. Firstly Newtonian gauge, unlike synchronous gauge, is well-defined and should yield unambiguous 
conclusions; and, moreover, the use of Newtonian gauge enables us to employ the easily extensible cmbeasy [59j code to 
track the size of the perturbations. We find that linear perturbations for a model Einstein de-Sitter universe introduce 
an effective energy density p e s ~ (4 x 10~ 5 )p m with an equation of state p e ff/p e s = w c s —1/19, while those for the 
ACDM concordance model introduce p G ff ~ (1-3 x 10 _5 )/? m with a marginally lower equation of state. Employing a 
modified Halofit [60( code to include quasilinear scales increases the effective energy densities to p e g « (5.6 x 10 _5 )p m 
and poff ~ (1.6 x 10 -5 )/5 m for the EdS and ACDM cases respectively, with unaltered equations of state. These are in 
good agreement with other evaluations of p c g (e.g. [H, HH H2] ) which typically remain of the same order. 

We begin with a brief overview of the 3+1 formalism of general relativity and in mil Al briefly discuss averages in 
cosmology before applying a simple average in ^IIIBI In $YV Al we calculate the forms of the backreaction terms for 
conformal Newtonian gauge and in £)V Al we present the deviations from an FLRW model from linear perturbations 
employing a modified Boltzmann code for both Einstein de-Sitter and ACDM universes, making it more general than 
the similar studies in [3l|, Hl[, who did not consider ACDM universes, and [HI, who worked in a rather narrow range 
of domain sizes. We follow by estimating the impact from quasilinear scales in gV Bl 



II. EINSTEIN EQUATIONS IN 3+1 FORM: BASIC QUANTITIES AND RESULTS 

We here briefly summarise the main results of the 3+1 formalism of general relativity and refer the readers to 
[U [111, 63] for details. Provided that the spacetime (A4, g) is globally hyperbolic we may foliate M. with a family 
of spacelike hypersurfaces S t , such that each hypersurface is a level surface of a smooth scalar field t, which we will 
later identify with the time coordinate. The normal unit vector n to the surface S t being timelike must satisfy 

nn= 1. (1) 

It is necessarily colinear to the metric dual Vi of the gradient 1-form di and the proportionality factor is called the 
lapse function a 

n = -aVt, (2) 

where a = (— Vi • Vt) -1 ^ 2 = — ((di, Vi)) 1 ^ 2 . To ensure that n is future oriented, a must be larger than zero and in 
particular never vanishes for a regular foliation. Since (di, an) = 1 the hypersurfaces are Lie dragged by the normal 
evolution vector an. 

Now introduce coordinate systems (x l ) on each hypersurface St that vary smoothly between neighbouring hyper- 
surfaces, such that (a; M ) = (i, x 1 , x 2 , x 3 ) form well-behaved coordinate systems on AL The 1-form di is dual to the 
time vector d t . Hence, just like the normal evolution vector, the time vector Lie drags the hypersurfaces. In general 
these vectors do not coincide. The difference defines the shift vector j3 

d t = an + (3. (3) 

It is purely spatial since (dt, (3) = (dt,dt) — (di, an) = 0. In terms of the natural basis of the introduced 
coordinates (x^) the normal unit vector can then be expressed as 

n" = ±(1,-13*). (4) 
a 

With the projection operator into St 



hfj,u = gp.v + n^riv (5) 



the components of the 3-metric hij become 

hij = g^hfti-. (6) 
From here one may find that the line element in these coordinates is 

ds 2 = (-a 2 + f3, t l3 l )dt 2 + 2/3 i dtdx i + h,,d.r'd.r J . (7) 



where /3, = hij/3- 
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The projection of the gradient V M of the unit normal vector defines the extrinsic curvature tensor 

Kii = -h%V k rij- (8) 



Ly = -n l v k nj 

Using Frobenius' theorem this can be expressed as the Lie derivative 



Kij — n^-nhij (9) 



and the evolution of the 3-metric hij is given by the Lie derivative along the normal evolution vector an. Using ([3|) 
and the evolution equation becomes 

aCnhtj = Cd t h i: j - Cphij = - 2T>^(3j) = -2aK i3 (10) 

where an overdot denotes d/dt, T>i is the covariant derivative on the 3-surface and brackets denote symmetrisation 
on the enclosed indices. 

Now we perform a 3+1 decomposition of the stress-energy tensor 

= pn^riv +2n( A J„) + Spa* p = T^n^rf , j t = -n M T iM , Sij = T%j (11) 

and project the Einstein equations onto the hypersurface and along its normal. This results in the Hamilton constraint 
equation 

K + K 2 - K)K{ = IGnGp + 2A, (12) 
where 1Z is the Ricci scalar on Et and K = K\\ the momentum constraint equation, 

Vj (K lj - h lj K) = BnGp; (13) 

and the evolution equation 

k i:j = a (Hij - 2KfK aj + KKij - %kGS %j + AirGhij (S - p) - Ahy) 

~ (ViVja - (3 a V a K l3 - 2K a(l V i) p a ) . (14) 

In line with much of the literature we set the shift vector to zero. 
The system now becomes 

ds 2 — —a 2 dt 2 + hijdx l dx^, 

-,0), n ll = (-a,0), K ij = -—h ij , 
a J la 

K + K 2 - K)K{ = WirGp + 2A, Dj (K lj - h ij K) = 8wGf, 

-Kij = Hij - 2K?K aj + KK l0 - SnGSij + ^Ghij (S - p) - Ah tJ - -V^a. 

III. AVERAGING IN COSMOLOGY 
A. The Averaging Problem and Averaging Procedures 



As we mentioned in section |H the Einstein equations have to be averaged over some domain to give an accurate 
description of the large scale behaviour of the universe. Since there is no physical meaning in the comparison of tensor 
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quantities at different spacetime points, appropriate averaging is very involved. Despite many promising approaches 
such as Zalaletdinov's macroscopic description of gravity [6J] and the regional averaging and scaling approach by 
Buchert and Carfora no procedure is yet agreed upon and work in this area is ongoing (e.g. [65l. l66l] ) . 
Lacking this process at current, the standard approach employs the following. Let 

V= I Vhd 3 x (15) 
Jv 

be the definition of the volume of a domain DCS smaller than the Hubble volume. Then one defines the average of 
a quantity A in this domain as 

{A) = 1 / AVhd 3 x. (16) 
v Jv 

This process has several drawbacks. Most importantly it is not invariant under a change of coordinates. Therefore, the 
averaged quantities lose their tensor character and moreover depend on the coordinate system in which the averaging 
was performed. Nevertheless, it is believed that this averaging process does make sense for scalar quantities. See 
[67J for a detailed discussion and an improved averaging process that produces quantities with approximate tensorial 
character in the weak field limit. To employ (|16p we will project all tensor quantities onto scalars. 



B. Average Hubble Rates, Scale Factors and Averaged Cosmological Equations 

We define an average Hubble rate by 

av = v = i d r ^ 3 = i r i ijh ^ 3 

av V VdtJ v VJ V 2 13 

Using the extrinsic curvature we thus have 

3— = - (aK) ; (17) 
av 

that is, the effective Hubble rate in a domain is given by the averaged trace of the extrinsic curvature in that domain, 
weighted by the lapse function. We can also quickly find the commutator between time and space derivatives, 

A ) = §1 {A) + V {A) + {AaK) = §i {A) + 3 S {A) + {AaK) • (18) 

Consider first the Hamiltonian constraint; we can immediately see that 

(a 2 TZ) + {{aK) 2 ^ - (aK)aK{\ = 16ttG (a 2 p) + 2A (a 2 ) . 
This is quickly converted into 



^ 2 = — (a 2 P ) + U* 2 )-UQv + Kv) (19) 



s av J 3 3 6 

with Buchert 's kinematical backreaction 

Qv = (a*(K*-KiKf))- 2 -{aKf (20) 

and 

TZv = (a 2 n) (21) 



a correction arising from the spatial curvature. Equation (|19p can be interpreted as an averaged Friedmann equation 
in the domain under consideration. 

The trace of the evolution equation for the extrinsic curvature is, 

aK - aKjaKj + aV l V l( x = AirGa 2 (p + S) - a 2 A (22) 
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which we have multiplied through by a 2 for simplicity of manipulation. Now, 

(aKj = {{olK))' + (a 2 K 2 ) - (aK) 2 - (aK) , 



and using the averaged Friedmann equation (|19p we find that the evolution equation averaged in the domain under 
consideration is 

— = ~ <« 2 (P + S)) + ± (a 2 ) + 1 -(Q V + V-o) (23) 

CLx> 6 6 6 

with 

Vv = (aDWiOt) - (aK) , (24) 

Buchert's dynamical backreaction. We interpret l|23p as an effective Raychaudhuri equation. Note that in the syn- 
chronous gauge, Vv = 0. 

Finally we have the so-called "integrability condition" which arises from the total continuity equation. We can find 
this by taking the time-derivative of the Friedmann equation and substituting in the Raychaudhuri and Friedmann 
equations. Writing Hx> = (H) where TL = — aK is a local Hubble rate, and using local matter continuity 

p + m(p+\s]=Q (25) 
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we can ultimately find the integrability condition 



d t {a^Qv) , d t (a 2 v (a 2 n)) I 2 f8nG A 



+167rG ((a 2 S) (H) - (a 2 SH)) + 6A {(a 2 H) - (a 2 ) (H)) - AH V V V - (26) 

It is easy to see that in dust- filled synchronous gauge models, with a = 1 and 5 = 0, the source on the right-hand side 
vanishes. For non-dust models, there is still an extra source term dependant on the local pressure and the local and 
averaged Hubble rates. In [451 ] . the authors make much use of the integrability condition, employing it iteratively to 
recover the backreaction at higher-orders in perturbation theory 1 . Since the authors employ synchronous gauge and 
we employ Newtonian gauge our results are not directly comparable, but theirs is an interesting approach and exploits 
the simple nature of the synchronous gauge integrability condition to extend the analysis as far as perturbation theory 
holds. 



IV. THE CONNECTION WITH LINEAR PERTURBATION THEORY AND THE SIZE OF THE 

BACKREACTION 



A. Backreaction Quantities in Newtonian Gauge 



The Buchert equations are exact for any inhomogeneous model; for further progress we must specify this model. An 
increasingly common approach is to employ a model that averages out to be Robertson- Walker on large scales but is 
a relatively realistic approximation to the local universe on smaller scales (for a recent study in this area sharing some 
of our aims see for example (HI). We choose instead to use a linearly-perturbed FLRW model, employing Newtonian 
gauge. While much of the previous study in this field has been undertaken in synchronous gauge (see for example 
[45j,|68|), this gauge contains unphysical modes and the metric perturbation can grow to be relatively large. While this 
is not necessarily an issue, Newtonian gauge is unambiguous and the variable <fi remains small and well-defined across 
almost all scales. Moreover, it is easily incorporated into Boltzmann codes and should serve as a complementary probe 
to the previous studies, although it should be remembered that we cannot directly compare calculations performed in 
different gauges. As we will not employ the inte grab ility condition, the complication in its source term is not for us an 
issue. It has been appreciated for some time [28l . |30| that the impact on large-scale evolution from perturbations is not 



1 See also |6Sl | in which the same authors employ the integrability condition to extend their analysis to scales of 50Mpc; this paper was 
completed at the same time as the first version of the current work. 
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expected to be large, something clearly seen in Newtonian gauge; since the perturbations themselves are consistently 
small across almost all scales the impact from backreaction will be at most of the order of 10 -5 . To our knowledge, 
however, this has not been fully quantified. It would be interesting in particular to study the behaviour of the 
corrections to the standard FLRW model as a function of redshift. 

We consider a universe filled with a cosmological constant and pressureless dust. This will serve as a reasonable 
approximation to the current universe on the largest scales. Working to first-order in the gravitational potential 0, 
we consider only scalar perturbations and work with the line-element 

ds 2 = -(1 + 2(f))dt 2 + a 2 (t)(l - 2(f>)S ij dx i dx j (27) 

where the potentials are equal as a pure dust fluid implies vanishing anisotropic stresses, and the scale factor a 
is distinct from the averaged scale factor ax>- We neglect the tensor perturbations as sub-dominant to the scalars 
although a full approach should naturally take these into account. When performing our averages across scales 
approaching the Hubble horizon, we consider the quadratic order in linear perturbations and neglect the averages of 
pure first- and second-order perturbations. Were we to consider averaging across relatively small scales then these 
neglected terms should also be taken into account for a full evaluation. Likewise we can neglect the second-order 
vector and tensor perturbations that scalars inevitably source, but again these should be taken into account for a full 
treatment on small scales, a problem complicated by the lack of usable covariant averaging procedures. Since we are 
considering purely scalar perturbations, all vector quantities can be written as the gradient of a scalar quantity; in 
particular, v % = —d tl ip for some velocity potential ip. 

Previous studies have occasionally extended their domain into super- horizon scales (e.g. [32l. l69l. l70l. FrTI. l72l. l73l l74j ) . 
This has been frequently criticised in the literature (e.g. [H, [7f| [76|, [z3] ) • While it is perhaps still arguable that super- 
horizon isocurvature modes might generate an observable impact (e.g. [Zil [zl, [zH ) ' we wm later restrict ourselves to 
adiabatic modes for which the effect is pure gauge and so choose our domain to be smaller than or equivalent to the 
horizon size. 

By unambiguously identifying the 3+1 and Newtonian gauge co-ordinates and expanding quantities up to the 
second-order in perturbations (see also [7^|), we can quickly see that 

a 2 = 1 + 20, a= V(l + 20) « 1 + 0- ^0 2 , o+ 1 «l-0+^0 2 , a~ 2 « 1 - 20 + 40 2 . (28) 

We also get 

hij = a 2 (t){l - 20)%, h ij «a- 2 (i)(l + 20 + 40 2 ),f+ (29) 
We can now proceed to evaluate the geometric quantities we need in the averaged equations. The extrinsic curvature 

is 

aK\ = - ^-0(1 + 20)") 5). (30) 

Using the Ricci scalar on the spatial hypersurface the curvature correction is 

Tlv = (a 2 Tl) = \ ((2V 2 + 3(V0) 2 + 120V 2 0)) (31) 

where V = (d/dx, d/dy, d/dz) as in standard vector calculus. 
From the extrinsic curvature we can immediately see that 

3^ = „( a X) =3 /^_0(l + 20)\ 
av \a I 

implying 

-- = £-(0(1 + 20)) (32) 

giving the relationship between the physical averaged scale factor and the scale factor employed in the perturbative 
approximation. If = then these coincide; this occurs in an Einstein de-Sitter universe, or when one considers a 
domain sufficiently small that its time variation can be neglected but sufficiently large that linear perturbation theory 
may be employed. (As an aside it should be commented in both these cases that feedback from the backreaction will 
rapidly break the condition = 0- even if this is at a negligible level.) 
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We can now rapidly find the kinematical backreaction, 



Qr = n( ( (^-#l + 2#M-/^-#l + 2#\ j =6((o-)-(„) ) CW) 



Using the covariant derivative on the 3-surface and 



aif = 30 2 - 3-^(1 -2(f>), (34) 



the dynamical backreaction is 



Vv = 3- (#1 - 20)) - 3 (</> 2 ) + ^ (V 2 ^ + 20V 2 - 2(V0) 2 ) . (35) 

This gives us the geometric information that we need. 

Consider now the fluids. We separated the stress-energy tensor with respect to the normal vector 

n"= (l-0+^ 2 ,o). (36) 

However, the stress-energy tensor employed in cosmological codes is defined with respect to the 4- velocity and for the 
case of pressureless dust is typically taken to be 

T,i V = eu^u v (37) 

where e is the energy density as defined in this reference frame, which will be linearised to e = p(l + 6) for an average 
p. The density p and FLRW density e thus do not coincide unless the velocity is orthogonal to the normal, which is 
not in general the case. The 4-velocity, to second-order in first-order perturbations, can be seen to be 

W+i(aV + 30 2 )x) , u»= (^-(^ + 4>+\{a 2 v 2 a 2 {l-2<t>)v)j (38) 
with Vi — Sijvi . Then the density is 



p = T^nV = e (u^) 2 =e(l + \a 2 v 2 ^j = e (l + aV) 



(39) 



Scaled by the lapse function, the quantity that enters the Buchert equations is then 

(a 2 p) = p (1 + S + 2cf) + aV + 24)5) (40) 

and so we define a "density correction" 

Tv = ^-p(S + 2<P + 2<f>5 + a 2 v 2 ). (41) 

The Friedmann equations can then be written as 

ax>\ 2 87rG A . „ ii-v 47rG A 

=— ^+3 +AF ' ^ + 3 +Ai? (42) 



with corrections 



AF = T v -l(Q v +K v ), &R=-\t v + \{Q-d + Vv)- (43) 



The effective pressure and density of the corrections are then 



^-Peff =Tv-\(Qv + K v ) , 16^Gp of f - \n v -On}- \vt> (44) 
3 3 3 
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and so the effective equation of state is 



1 (TVd- 4V v -3Q-d \ 
W ° S = -A n v -6T v + Ql '■ (45) 



Note that these dehnitions differ from those presented before due to the presence of the T-p term, arising from directly 
comparing with a reference FLRW model. 

We take the domain to be large enough to allow us to neglect the averages of first order quantities on the background. 
Specifically, we take the domain size to be of the order of the comoving Hubble scale, across which first-order averages 
can consistently be neglected. While this implies that our analysis can only be taken to apply on the very largest 
scales, taking the domain to be so large also implies that we can invoke the ergodic principle to convert our spatial 
averages into ensemble averages; we exploit this in the next section. On such large scales, the modifications become 

n v = 4<(V0) 2 + , (46) 

' 12 



Vv = 6(F), (47) 

Q v = (#) - 3 (0 2 ) + 1 {4>V 2 <t> - (V</>) 2 ) , (48) 

Tv = ^P(2#+«V). (49) 

We can see that, depending on the signs of (4>V 2 (f)) and (4>8), for an Einstein de-Sitter universe with = (as 
considered in 43]), one might get either an enhancement or reduction of both the effective Hubble rate and of the 
effective acceleration. 



B. Ergodic Averaging 

The simplest attack on the above formalism is to employ a Boltzmann code such as cmbeasy or cmbfast [80j . 
The use of linear theory automatically implies that we are working on relatively large scales; we can render the 
system developed in the last section more tractable by taking our domain to approach the Hubble volume itself - 
large enough, at least, that we can employ the ergodic theorem and convert the spatial averages into averages across 
a statistical ensemble. (Similar approaches were also employed in [28|, |30|.) 

Consider first the general average (A(t, x)B(t, x)) where {A, B} £ {<fi, 5,v, cb}. Then 

(A(t, x )B(t, x )) = (A(t, k)B* (t, kO) e-« k - k >■* . (50) 

At linear order and in the absence of decoherent sources, different wavemodes are decoupled from one another and 
the evolution equations depend only on the magnitude of k and so we can write 

A(t,k)=a(k)A(t,k) (51) 

for any quantity. The initial correlation between two arbitrary quantities is 

(«(k)U'(k')) = ^(#-k') (52) 

where V^,{h) is the primordial power spectrum of the metric fluctuations. Inserting this into the correlation and 
integrating once yields 

f Ah 

(A(t,x)B*(t,x))w J V^k)A{t,k)B*{t,k) — . (53) 
The averages involving gradients are only a little more complicated; for ((V<^) 2 ) we have 

<(V^) 2 ) - | k | k/ Hk)(zkO(0(t,k)0*(t,kO> e ^ k - k > x ^^ (54) 

= f k 2 r^(k)m,k)\ 2 j-, (55) 
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which is positive-definite as it should be. The other curvature term is 

= - J ev t {k) i«t,fc)i 2 *. 

This term is then negative definite. 

The kincmatical, dynamical, curvature and density corrections (|46H49[) . then become 



(56) 
(57) 



Qv 



6 / V^k) 4>{t,k) 



2 dk 



= -3 



8ttG 



V(fc) ( j^t, fc)f + ~ fat, k)<j>*(t, k) + f(t, k)j>(t, kj) + ~ m, k)\ 2 \ ^ 



dk 



i> \ I k 2 V^{k)a 2 {t)\v{t,k)\ 2< ^~ I V^{k){ct>{t,k)6*{t,k) + (j>*{t,k)8{t,k))^- 



k 



dk 
~k 



(58) 
(59) 
(60) 
(61) 



For a more realistic approach where we separate the baryons from the cold dark matter one can readily see that the 
contribution from the density can be expanded as 



8nG 



k 2 r^(k) (p b \v b (t,k)\ 2 + p c \vc(t,k)\- 



dk 



dk 

V^k) m, k) (p b S* b (t, k) + p c 5* c (t, k)) + <j>*(t, k) (p b 5 b (t, k) + ~p c 5 c {t, k))) —. 



(62) 



Since Boltzmann codes tend to be written in conformal time we must these into conformal time expressions using 
d/dt = (l/a)d/dr)} denoting conformal time quantities with a tilde we thus have 



4>' a _ 1 a 1 
a ' a a a 



(63) 



in the expressions above. 



V. RESULTS 



A. Corrections to the FLRW Picture from Linear Perturbations 



We incorporate the above formalism into the cmbeasy Boltzmann code and run consistency checks with the cmbfast 
code, modified to output conformal Newtonian quantities. We will evaluate the terms across all post-recombination 
redshifts and imploy and infra-red cut-off at the comoving Hubble scale, fc m i n = 1/r), avoiding the unphysical gauge- 
dependent super-horizon contributions. The small-scale limit fc max of our domain, determined from the stability of 
the integration with respect to changing fc max is about fc max ps 30MPc _1 and we integrate to fc max = lOOMpc -1 . 
Naturally, we do not claim that our results at such small scales are complete, merely that we are evaluating the 
contribution of such large-scale modes on these scales. 

At a linear level, there are scaling relations that hold to a high level of accuracy; it is obvious from Poisson's 
equation that 



(f) oc S/k 2 ; 

from Euler's equation v + (d/a)v oc V0 we can also predict that 

\v\ oc S/k. 



(64) 



(65) 



From here, one may immediately state that on smaller scales where 1/fc 4 <C 1/fc 2 and for an approximately scale- 
invariant primordial power spectrum, 



S 2 S 2 
Qp(fe) oc -tj, Tt>(k) oc Vv{k) oc TZv(k) oc 



(66) 
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where Qp(fc) is the integrand of Q-p. Qx> is thus generally subdominant to the other corrections except on very large 
scales. This also implies that each correction is approximately of the form 



A D =aJ \S\ 2 ^dk (67) 



for some constant a, reminiscent of the approximation in equation (32) of [28| . 

For our models we take a low-Hubble constant Einstein de-Sitter model with f2& = 0.05, fi m = 1 and h = 0.41 
and a WMAPIII ACDM concordance model. The left panel in Figure [T] shows the (integrands of the) four different 
correction terms at redshifts of z = 10 and z — for the EdS case and the right panel the same for the ACDM case, 
and in Figure[2]we present the correction terms as a function of z for the EdS and ACDM cases; the premultiplication 
by a 3 (z) acts as a volume normalisation. The subdominance of Qx>(k) and proportionality between 1Zt>, 7v and "Pp 
is very clear, with TZp the strongest correction. 

From the proportionality of the other corrections (which holds up to relatively large scales) we can write Tp = tlZz> 
and Vp = pTlv and express the equation of state (|4"5)) as 

«*«4(^V (68) 



3 V 1 - 6t , 

which is simple to evaluate numerically by selecting some pivot scale on which to evaluate the ratios p(z) and t(z). 
Doing so at the current epoch and a pivot scale of k = O.OlMpc -1 , we find p w 2/9 and t « 1/20; both of these are 
also only slowly evolving. This gives an estimate of the average equation of state as w e g « —1/19. In the left panel 
of Figure [3] we plot the evolution of w e s for both the model EdS and ACDM cases evaluated directly from equation 
fl4"S")) . In both cases the effective equation of state from linear perturbations remains around the order w e s w —1/19 
and the correction terms to the usual FLRW as a whole thus act as a form of non-standard dark matter and not as a 
dark energy. These results compare reasonably well with the estimate in [28[ where the author found w e g w — 1/27 2 . 
To quantify the impact on the Friedmann and Raychaudhuri equations, define 

_ 8ttGp„ i0 _ AirGp~ m0 

as the standard contribution from baryons and CDM. We can then consider AF/F m and AR/R m as a well-defined 
measure of the impact of the correction terms in both EdS and ACDM models. By equation (|44]) . AF/F m for an 
EdS model is just p c g/'p m . For the sake of clarity we choose to focus on deviations from the standard behaviour but, 
naturally, we can re-express these quantities in terms of the effective energy density and pressure. 

For the EdS case we can see in Figure [3] that the total impact from linear perturbations tends to ~ 4 x 10~ 5 on 
the Friedmann equation and ~ —3.2 x 10~ 5 on the Raychaudhuri equation; the corrections thus act with a positive 
effective density and with insufficiently negative pressure to accelerate the universe and instead act to decelerate it 
at a negligible level. The dashed curves in Figure [3] are for the ACDM case; they tend to impacts of ~ 1.3 x 10~ 5 on 
the Friedmann equation and ~ —1 x 10~ 5 on the Raychaudhuri equation. The behaviour, as might be expected, is 
qualitatively similar to that in the EdS case and the impact is significantly less (by a factor of roughly 10/3 at the 
current epoch). 

This compares reasonably well with the approximation in [28| which in our notation is 

^ = ^VGPJ / ^dk. (70) 

Pm 4 J k 

Evaluating this approximation for the EdS universe, we find that this [under/over] estimates the correction terms with 
respect to our more detailed study by a factor of about 30%. In the same paper, the author found a "cosmic virial 
theorem" , 

Pm = "-Poff (71) 

where p m is the correction to the matter pressure arising from gravitational interactions. Crudely modelling the 
matter pressure as p m » (l/3)p TO (w 2 ) we find that p m underestimates — p e g, again by ~30%, but the proportionality 
between them holds remarkably well for a wide range of redshifts. 



They also found w e g ~ —1/15 for a clustering cosmon field but we have not considered such a component here. 
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FIG. 1: The correction terms T-p, Vt>, TZ-t> and Qt> as a function of k at z = 10 (black) and z = (grey) for (left) the sample 
Einstein de-Sitter model and (right) the WMAPIII concordance model. 
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FIG. 2: The evolution of the corrections for (left) the sample EdS model and (right) the WMAPIII concordance model. 



For the ACDM case we can consider an alternative normalisation that directly quantifies the impact on a universe 
with a cosmological constant, 

AF AF AR _ AR 

d-o/a-D F m +A/3' \av/av>\ \Rm ~ A/3| 

The benefit of doing so is that it demonstrates the declining contribution of the corrections with respect to the cosmo- 
logical constant; however, it also introduces a singularity at a redshift of z « 0.8 when the reference FLRW universe 
undergoes a transition from deceleration to acceleration. The impact on both the Friedmann and Raychaudhuri 
equations tends to ~ 10~ 6 and the maximum contribution is at a redshift of z « 1.3. 
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FIG. 3: Left: The effective equation of state w e g as a function of redshift for a sample Einstein de-Sitter and WMAP concordance 
models. Right: The impact of the corrections onto the Friedmann equation (black) and the Raychaudhuri equation (brown), 
normalised in the ACDM case to the matter content (dashed) and to the standard equations (solid). 



B. Corrections to the FLRW Picture from Quasilinear Perturbations 

In the previous section we evaluated the impact on the large-scale evolution of the universe from linear perturbations, 
demonstrating that it is small, as should be expected, but also that it acts as a dark matter and not a dark energy. 
In this section we employ the halo model to estimate the impact of perturbations on smaller, quasilinear scales which 
are much less understood. 

The publicly-available Halofit code [6(| converts a linear CDM power spectrum into the CDM power spectrum 
from the Halo model. We employ a modified Halofit code that instead takes the linear matter power spectrum and 
estimates the nonlinear matter power spectrum. The square-root of this power spectrum can then be used to estimate 
S(k) and we can then recover the other relevant quantities from the scaling relations v oc S/k, <f> oc S/k 2 . If we define 

/(*) = Pm (^) 

we can employ this as a scaling factor to recover estimates for the quasilinear behaviour of 5(h), v(k), <j>(k) and 
4>{k). We can then employ the formalism we developed for the linear case to estimate the impact in the quasilinear 
case, retaining the same domain size. We should immediately note two caveats. Firstly, while the scaling relation 
for the velocity is extremely precise for linear perturbations it is not for non-linear perturbations and so we have 
automatically introduced a source of error. Perhaps more importantly, we are employing a formalism developed for 
linear perturbations in which different wavemodes decouple from one another, allowing us to separate the statistical 
problem into transfer functions and a primordial power spectrum. This does not hold for a nonlinear problem. For 
both of these reasons, our results are only intended to be taken as good approximations on quasilinear scales until 
the velocity virialises, at which point a more detailed study is necessary. 

In Figure 2] we present approximations for the corrections to the ACDM model at redshift 0. As before, the 
kinematical backreaction is strictly negligible and as we have merely scaled the previous corrections by the same 
quantity, the effective equation of state remains unchanged. The total impacts on the Friedmann and Raychaudhuri 
equations for both are shown in Figure 0J We see that although the halo model provides a boost in power on smaller 
scales, and even though TZx> and Tt> in particular include factors of k 2 in the integral that increase the contribution 
from smaller scales, the impact is not significantly greater than from linear scales. More quantitively, for an EdS- 
universe the impact at the current epoch on the Friedmann equation is of order 5.6 x 10 -5 and on the Raychaudhuri 
equation is of order —4.7 x 10~ 5 ; normalised to the matter content, the impact at the current epoch from ACDM 
perturbations is of order 1.6x 1CP 5 on the Friedmann equation and of order — 1.4x 1CP 5 on the Raychaudhuri equation, 
normalised to the full evolution 4.4 x 10~ 6 and —2.7 x 10~ 6 , respectively. 
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FIG. 4: Left: The halo model approximations to the quasilinear corrections for the WMAPIII concordance model at z = 0. 
Right: The halo model approximations to the quasilinear corrections to the Friedmann (positive) and Raychaudhuri (negative) 
equations; the linear prediction is in grey. 



VI. DISCUSSION AND CONCLUSIONS 



We have presented quantitive estimates of the corrections in Newtonian gauge to a standard perturbative FLRW 
model from an explicit averaging procedure, which can be separated into four distinct quantities: the kinematic 
backreaction, which remains strictly negligible across the scales of interest, the dynamic backreaction, the impact of 
spatial curvature and a correction to the FLRW energy density. In line with expectation, the impact from linear 
perturbations is insignificant and of the order of 10 -5 for both ACDM and an Einstein de-Sitter model. Specifically, 
the effective energy density arising at the background level from inhomogeneities is p e g ~ (4 x 10 _5 )p m and p c g ~ 
(1.3 x 10 _5 )p m for the EdS and ACDM models respectively. Moreover, we have presented estimates arising from halo 
model corrections on quasilinear scales which are not much larger than those from the linear perturbations, the largest 
contribution arising from the quasilinear modes in an EdS universe remaining below 10 -4 ; we find p c s « (5.6 x 10 _5 )p m 
for the EdS model and p c g w (1.6 x 10 -5 )/J m for ACDM. This is in broad agreement with recent calculations from 
Vanderveld et. al. estimating the total possible backreaction from observations of Type la supernovae and agrees with 
a growing consensus that the impact of backreaction across all scales remains at the level of 10~ 5 — 10 -4 . Additionally, 
the total effective equation of state arising from the different modifications does not act as a dark energy and instead 
as a dark matter with a slowly-varying equation of state w e g « —1/19 for z £ (0, 10). The impact is to decelerate 
the universe and, contrary to other studies, the corrections impact positively on the Friedmann equation implying a 
positive effective pressure. (The discrepancy arises due to our definitions of the effective energy density and pressure 
of the corrections taking into account the differences in definition of energy density in an FLRW frame and with 
respect to our foliation.) 

This does not, however, necessarily imply that there are no observational consequences arising from inhomogeneities, 
particularly on small scales, as our analysis is limited to the impacts in very large volumes and breaks down when 
mode-coupling and virialisation become significant. In particular, one can readily imagine situations in which the 
local Hubble rate is significantly larger than the global average (the so-called "Hubble bubble"). The impact of 
local inhomogeneities directly on luminosity distances has also long been well studied and remains an active area 
of interest. Moreover, the issue of averaging in cosmology and general relativity in particular remains very much 
an active field. The average we have employed in this study is certainly not perfect; even were we to extend our 
approach to second-order perturbations it cannot cope with the vector and tensor perturbations that necessarily 
become significant. Improving the average and applying it to nonlinear situations is an immediate aim. 

The prospect for future study remains large. Our model here is necessarily simple and there are many avenues both 
analytical and numerical that can be pursued. Much interest is currently focusing on the exact swiss-cheese models 
and particularly on the impact on direct observables from these models. The prospect of directly calculating both 
the direct backreaction in different domain sizes and the impact on null geodesies (and hence luminosity distances) in 
different matter configurations from numerical simulations is also an exciting one. More immediately, one can imagine 
employing second-order perturbation theory in which the ambiguity between averages taken across the assumed zeroth- 
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order manifold and those across the linearly-perturbed manifold may become significant; depending on the surface 
one defines a vanishing (</>(2)) across, one would expect a non- vanishing contribution on the other. In particular, if one 
defined second-order perturbations to vanish on average on the background then there would be a residual contribution 
to the linear perturbations even before any other backreaction is taken into account. This is an interesting issue that 
deserves more study. More ambitiously, the fully non-linear approaches of Vernizzi and Langlois (8ll . 182 . |83| and 
Enqvist et. al. [84TJ might provide an interesting tool for studying the impacts of averaged nonlinear perturbations. 
The direct backreaction from more complicated systems than our linearly-perturbed FLRW universe might also be 
considerable; an immediate example, first considered in [28| . would be the backreaction from inhomogeneities in a 
cosmological scalar field which were shown to potentially have an impact of order unity. Such systems might also be 
considered in a fully non-linear manner [85| . 

In summary, one may comment that backreaction and the impact of inhomogeneities are direct physical effects and 
should certainly be taken into account in any full treatment of the observables. However, the direct impact when 
evaluated for the concordance cosmology on very large scales is too small to provide an alternative to the dark energy. 
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